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We suggest the possibility of creation in the early Universe of stable domains of 
radius a few kilometers wide, formed by coherently excited states of 7r-mesons. Such 

x: 

domains appear dark to an external observer, since the decay rate of the said coherent 



pionic states into photons is vanishingly small. The related thermal insulation of the 
domains from the outer world could have allowed them to survive till present days. 
The estimated maximum radius and the period of rotation of such objects turn out 
to be compatible with those of certain pulsars. 



00 i I. INTRODUCTION 

O 

Low-energy QCD can be characterized by four nonperturbative quantities, of which the 
gluon condensate (g 2 (F£ v ) 2 ) = (G 2 ) and the vacuum correlation length a are related to 
confinement (so that the string tension between the two static quarks in the fundamental 
representation is Of oc a 2 • (G 2 ) while the chiral condensate (ipip') and the constituent 
quark mass m characterize spontaneous breaking of chiral symmetry. One can still consider 
an extrapolation of the chiral condensate (z/>V) to the heavy-quark limit (where (i J ' l l J ) he£lvy 
ceases to be an order parameter of spontaneous chiral symmetry breaking). In this limit, the 
role of the constituent mass of the quark is played by its current mass M, yielding the known 
formula of QCD sum rules j^j]: {i J i J ) hea:vy ^ — (G 2 ) /M. In fact, as long as the mean size of 
the quark trajectory in Euclidean space remains smaller than a, the nonperturbative Yang- 
Mills fields inside the trajectory can be treated as constant, leading to the area-squared law 



for the quark's Wilson loop 



l|. As has been shown in Ref. J3], this law reproduces correctly 
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the above-mentioned heavy-quark condensate \'0'0) lieavy > With the decrease of the current 
quark mass, the mean size of the quark trajectory gradually gets larger than a, and the 
area-squared law of the Wilson loop gets morphed into an area law. The effective scale- 
dependent string tension entering the area law characterizes the strength of the confining 
force between the light quarks much in the same way as the constant string tension does for 
heavy quarks. 

Alternatively, one can use Nambu-Jona-Lasinio-type models with confining kernels, 
henceforth called Generalized NJL (GNJL) models, as a tool to study the interplay between 
confinement and spontaneous chiral symmetry breaking. Unlike the standard Nambu-Jona- 
Lasinio (NJL) models [4| with either local or short-ranged kernels, which do not describe 
confinement, the latter is present explicitly in GNJL models [s - T ] , which successfully ful- 
fill the known low-energy theorems of Gell-Mann, Oakes and Renner Goldberger and 
Treiman [9|, Weinberg 10], and reproduce the so-called Adler zero 11]. Yet another positive 
property of such models is the stability of the said results against variations of the parame- 
ters, as well as the functional form of the confining kernel. Therefore we are led to conclude 
that these models do provide a robust mechanism for the existence of a Goldstone boson 
- the pion — which becomes massless in the chiral limit. The presence of this boson can 
be associated with the specific realization of spontaneous chiral symmetry breaking, which 
is related, through the confining kernels, to confinement. This mechanism is enshrined in 
what is known as the mass-gap equation (MGE). Solving this equation amounts to a proper 
choice of the quark Fock space J 7 , that is, to a proper choice of the fermionic ground state 
of the theory. In GNJL models, the MGE turns out to possess more than one solution, the 
replicas [6]. Among those, the one with the lowest energy yields the QCD chiral condensate. 
Higher replicas describe excited states, in which the absolute value of the chiral condensate 
is smaller (that is, the chiral symmetry is "less broken") than in the QCD vacuum. There- 
fore, one can naturally pose a question as how and where such excited states could exist in 
Nature. Because theoretical evidence for the possible existence of such states (which, as we 
shall see, can be realized as coherent states of pions) was restricted to GNJL calculations, 
it is important to verify, within the effective-action approach, which is a gauge-invariant 
formalism, whether such states may remain a degree of freedom of QCD. To this end, we 
study in Sec. II the correspondence between the effective-action approach of Ref. [3] and 
GNJL models. In this paper, we show that, despite their different nature, both approaches 
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turn out to possess, in the chiral limit, a scale-dependent effective mass, which is necessary 
condition for the existence of replicas. In the heavy-quark limit, the two approaches have 
already been shown equivalent, cf. Ref. Jsj]. 

As for the actual problem of possible realization of replicas in Nature, we argue that such 
replicas, once thermally excited in the early Universe, may form certain stable macroscopic 
domains. To explore such domains quantitatively, we calculate in Sec. Ill their internal en- 
ergy. It is given by the difference between the internal energies of the replica-filled domain 
and the physical vacuum outside that domain. As such, this internal energy is positive- 
definite, making the domains stable against the gravitational collapse. Owing to this fact, 
the domains can be studied by using the well-known Tolman-Oppenheimer-Volkoff equa- 
tion [12 1 . For sufficiently low temperatures, we show the pionic contributions to the replicas' 
internal energy to be negligible when compared with the "excited vacuum" contributions, 
thence allowing us to assume the internal energy inside such domains to be constant. 

This assumption allows us to apply the known form for the metric inside spherically- 
symmetric objects with a constant energy density (as provided by the corresponding solution 
to the Tolman-Oppenheimer-Volkoff equation 12|) to calculate quantitatively the leading 
gravitational correction to the chiral condensate inside the domains. As a result, we find 
that this correction could only become relevant at the distances < 0.55 fm to the center 
of the replica-filled domain, which are, for the first excited replica, 20 orders of magnitude 
smaller than the size of the domain. Such irrelevance of the gravitational correction to the 
chiral condensate makes our approach self consistent, since it justifies our use of the flat- 
space internal energy of the first replica. Furthermore, by a simple estimate we show that 
the anomalous decay-rate of coherently excited 7r-mesons, which fill the domain, into two 
photons is vanishingly small. Using this fact, we put forward a conjecture that such domains 
in the Universe have had a chance to survive till the present time. 

The paper is organized as follows. Section II is devoted to a theoretical discussion of 
replicas. We show how a long-range confining kernel can be modeled, in both GNJL and the 
effective-action approaches, by a scale-dependent constituent quark mass, which, in turn, 
becomes the cornerstone for the possible existence of replicas. In Sec. Ill, we calculate 
the internal energy of the domain filled with coherently excited 7r-mesons, as well as its 
maximum radius compatible with the stability against the gravitational collapse, and a 
possible period of rotation. In Sec. IV, we use the effective-action approach to evaluate 
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gravitational corrections to the chiral condensate, and show such corrections to be negligible. 
In Sec. V, we summarize the main results of the paper. In Appendix A, we recall the world- 
line representation of the Abelian Higgs model with a nonlocal quartic Higgs interaction, 
that is the bosonic prototype of an NJL model with a nonlocal kernel. In Appendix B, 
for the sake of completeness, we collect well-known results on the gravitational field of a 
spherically-symmetric object with a constant energy density. 



II. REPLICAS 



The possible existence, in QCD, of new types of physical states, hitherto named replicas, 
has been put forward by phenomenological studies based on the GNJL models ^, [y, ["LsL 



151 ] . In this paper, we follow two pronged approaches to the issue of replicas. First, in 



Subsec. Ill Al below, we give a bird's eye view on the physics and the origin of replicas 
within GNJL models. Then, in Subsec. Ill Bl we proceed to establish the connection between 
the GNJL and the effective-action approaches. This is important not only because of the 
relevance of validating, within a different setting, the possible existence of replicas, but also 
because, once such states are to fill a macroscopic domain (and hence to be a source of 
gravity), it is of paramount importance to study the back- reaction, which the gravity might 
have on the quark condensate (that has produced, in first instance, such gravity forces). It 
happens that this study can reliably be carried out within the effective-action formalism - 
see Sec. IV - which is not the case for GNJL models. 



A. The microscopic analysis of spontaneous chiral symmetry breaking and the 

origin of replicas 

In order to illustrate the correspondence between the GNJL and the effective-action 
approaches to chiral symmetry breaking, let us start with a brief introduction to GNJL 
models of QCD. Upon integration in the QCD partition function over the gluonic fields in 
the Gaussian approximation, that is, disregarding the explicit contributions of all connected 
averages (q n A°; 1 Af? ■ ■ ■ Af") with n > 2, one obtains an effective action containing a 
four-quark interaction: 



S = T 



<K 7 ^ + m)z/> + ~ / jZ(x)(g 2 A°(x)Al(y)) ( .,;:<,/) 



(1) 
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Here T stands for time ordering, = ip , y l jT a ijj is the quark current, T a is an SU(iV c )- 
generator in the fundamental representation, and we use a shorthand notation J = J d^x. 
Of course, the heavier the current quark masses m are, the better the Gaussian approxima- 
tion becomes, in which case the two-point connected correlation function {g 2 Af 1 (x)A b (y)) conn 
becomes the usual gluon propagator (g 2 Af l (x)A b (y)). For lower current quark masses, 
{g 2 Af l (x)A b / (y)) must be thought as representing an effective, phenomenologically given, 
force between quarks. Let us denote 

Df u (x,y) = (g 2 A^x)A b u (y)), m(x) = m + 2 [ D%(x, y)^T a (1>{y)${x)) m 1 V T\ (2) 

Jy 

where (. . .) stands for the Wick contraction. Using the Wick theorem, we can then cast 
Eq. (pQ) to the form 

S = S +: f$ M + m(x)} ^ + jZ(x)D; b u (x, y) 3 b {y)-- , (3) 



where Sq is a c-number resulting from the Wick contractions, and : : stands for normal 



ordering. There exists a vast literature devoted to phenomenological ansatze for D^(x, y) 



15, 



6|, aimed at the simulation of confining forces between quarks (and hence at the simulation 
of the contributions of the entirety of the cumulant expansion to the quark confining forces). 
It turns out that, regardless of the form of Dff v (x,y), we have the effective quark mass 
rh(x) ~ M for the heavy-quark current mass M. However, when the current quark mass 
m — > 0, fh(x) does not in general go to zero but rather becomes a function of the quark- 
antiquark separation. GNJL models provide one with a natural tool for studies of this 
phenomenon. Indeed, considering an instantaneous correlation function D^ b u (x,y) given by 
a suitable confining potential, one can apply to the problem the whole variety of Hamiltonian- 
based methods known in the literature. 

The second-quantized Hamiltonian corresponding to action (TjQ) contains anomalous (also 
called off-diagonal) terms. Among those, the bilinear anomalous terms coming from 
■ f $ bV^M + ip'- are of particular importance. In order to eliminate such terms, one 

uses a Bogoliubov-Valatin transformation V<B[<£> P ] [6] to represent the original quark fields 
via new "dressed" quark fields. The transformation Vi3[v? p ] can be thought of as a rotation 
in the Hilbert space spanned by the u and v spinors (which together compose a bispinor 
if}) parametrized by a momentum- dependent angle ip p . Within GNJL models, one further 
eliminates the anomalous terms by means of the mass-gap equation (MGE[<p p ] = 0), whose 
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solution(s) define the chiral angle ip p . It is therefore quite clear that the prerequisite for hav- 
ing non-trivial solutions for (p p is the existence of a non-zero m(x). The remaining quartic 
anomalous terms can also be eliminated by another Bogoliubov-Valatin transformation - 
see the discussion at the end of this subsection. 

Now, the mass-gap equation, being nonlinear, may have multiple multi-node solutions 
for the chiral angle (f p - the replicas. For the quadratic kernel, the replicas were discussed 
already in the pioneering papers {(]], while the long-range kernels of a generic form and 2D 



HQ and Ref . [15f , 



respec- 



QCD (the 't Hooft model) in this context were studied in Refs. |13|, 
tively. Existence of multiple solutions to the MGE was also confirmed by the calculations 
done in other approaches 16]. The zeroth replica y? P i?0 ' ) provides the vacuum state |0) with 
the lowest energy-density, and defines the corresponding fermionic Fock space J-q. The vac- 
uum energy-density corresponding to any other Fock space T n is larger. For J^q, the MGE 
coincides with the pion Salpeter equation and ensures, at the same time, both the chiral 
pion to be massless and n — n scattering lengths to be zero. In other words, it is in general 
unavoidable to have a chiral pion plus all the low-energy theorems thereto, and not to have 
an n-numerable set of solutions f^(p p Rn \ (vpip)^ to the MGE. 

An operator creating a certain replica R can be written in a closed form as 14 1 



\R) = S\0), 5=JI 

c,/,P 



cos 2 ^ + sin^cos^t p + lsin^ 
2 2 2 cJ ' p 2 



Aif p 



2 °c,/, P 



(4) 



where |0) and \R) are the physical BCS vacuum and the vacuum of the replica, which 
are defined by the chiral angles y?p°' ) and (fp , respectively, and Aip p = — (p p \ The 
operator & * creates a 3 P quark-antiquark pair of color c, flavor /, and with the relative 
3-momentum p: 

4c/ (p) 



C c,f,p 



[^/(p)^i/(p)]^3 Po 



M: 



Po 



{o-p)io- 2 , 



(5) 



where cr's are the Pauli matrices. It turns out that, for a sufficiently large spatial volume 
V^ 3 \ the states |0) and \R) are essentially orthogonal to each other [jjj]. That is because 
their overlap gets exponentially suppressed with the increase of the volume V^ 3 \ namely 



(0\R) = exp 



N r N f V {3) 



d 3 p 
(2^)3 



In ( cos' 



A(f p 



(6) 



The same is obviously true for the overlap of any two quark states belonging to different 
Fock spaces. 
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Furthermore, as was shown in Ref. 



131 ]. for any chiral angle respecting the MGE, yet 



another Bogoliubov-Valatin-like transformation can be performed in order to reformulate 
the GNJL model in terms of the observable mesonic states. Such a transformation can 
be applied to replicas as well, provided the volume is large enough. As a result, the 
quark BCS states \q)o built on top of the vacuum |0) get transformed to the mesonic states 
|M) built on top of the mesonic vacuum \Q)o, while the replica quark states \q)ji get 
transformed to some coherent mesonic states \M)r — excitations over \Q)r. It turns out 
that the operator S defined in Eq. (j3J) can be represented in terms of the compound mesonic 
operators. Since the dominant role in the cloud of coherent mesonic excitations is played 
by the lightest hadrons, it is natural to approximately consider mesonic replica-states as 
being formed by "almost non-interacting" pions. The weakness of pionic interactions is a 
consequence of chiral symmetry. Namely, as was already mentioned, for massless quarks, the 
7i — 7i scattering length vanishes, and the pions at low momenta appear almost transparent 
to each other 7j. With this remark, we finish the physical discussion of replicas in the 
context of GNJL models. In the next subsection, in order to evaluate the scale-dependent 
constituent quark mass rh, we resume this discussion within the effective-action approach. 



B. A correspondence between the GNJL and the effective-action approaches 

From now on, we work in Euclidean space and denote rh as just m, while reserving the 



symbol M for the heavy-quark current mass. In Ref js], we established the correspondence 
between the GNJL and effective-action approaches. This correspondence becomes exact in 
the heavy-quark limit, which is the limit where both the Gaussian and the minimal-area 
(squared) approximations become exact. For lower current quark masses, both approaches 
have to be supplemented by phenomenological-kernel ansatze. For the effective-action for- 
malism this necessity is clear: When we depart from the heavy-quark limit, that is, from 
small Wilson loops of sizes of the order of the correlation length a, we depart from the 
minimal area they span. This leads us from the minimal plaquette that supports one homo- 
geneous flux to an ensemble of correlated plaquettes, each having its own constant flux, which 
needs not to be the same across the neighborhood. Therefore, to establish a correspondence 
between the GNJL and the effective-action approaches, away from the heavy-quark limit, 
where they are proven equivalent, is the same as to establish a correspondence between their 
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respective phenomenological extensions for lower quark masses. To this end, and as a first 
step, let us revisit the heavy-quark limit, where the sizes of quark trajectories, being of the 
order of the inverse current quark mass 1/M, lie inside the circle of the size of the vacuum 
correlation length a. At distances smaller than a, the non-Abelian field-strength tensor F£ v 
can be treated as constant, yielding the gluon propagator 



K( x > y) * zwt> {fK, *%, {-^ ))■ (7) 

Furthermore, the correlation function of gluonic field strengths in this formula reads 

9 2 K, (^) F% (^) ) = (G 2 ) ■ j^j ■ ^{8^ - <W<W)- (8) 

Inserting Eqs. (j7]) and ([8]) into Eq. ([3]) and introducing an antisymmetric tensor current 
J® = \{jtx v — jfafj), we obtain the following GNJL action: 

s = JjM + M)4> + 96( ^1 1} I JU X ) J Uv) ■ ^(M^ - 5 -<W)- ( 9 ) 

This model has a confining kernel, which allows for replicas Notice that, in Ref. an 
extrapolation of the model to the chiral limit has been addressed. 

Now, instead of solving the MGE corresponding to the GNJL model, Eq. ([9]), one can 
map this model onto the one-loop quark effective action jjj. To do so, one notices that 
the four-quark interaction in Eq. ([9]) can be bosonized by introducing an auxiliary Abelian- 
like field A®(x) = \x v n a T v ^ where T vll is a constant Abelian strength tensor, and n a is a 
constant unit color vector, n a n b = 5 ab . This observation allows us to rewrite e~ s as 

e S = ) ( 10 ) 



where 



+oo 



and the covariant derivative corresponding to the auxiliary field A't reads = + 

iA a ^T a . The heavy-quark condensate is given by the tadpole diagram as 

' heavy = ~ \ tr iMA] + M 
It can be obtained upon the differentiation over M of the one-loop effective action of a heavy 
quark in the auxiliary field A a ^. 

<^> heavy = (tr In (iMA] + M)) T . (13) 



(^)hea vy = - < tr ^T7V-TJ ) ■ (12) 
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This expression is defined up to an overall normalization constant, which can be fixed upon 
a comparison with the free case, where the quark condensate vanishes. To divide the free 
case out means to perform the following substitution: 

tr In (j,D,[A] + M) -+ tr In 7 ^j^ M (14) 

We can now use the fact that the effective action of a quark in the field is known in 
an analytic form (see, for example, the two reviews quoted in Ref. }18|). This observation 
yields 

'Y D \A]-\-M f°° e ~ M2fi 
tr ln V^i^J ^ = _ 2N{ . tr r T a T a^ . / ds | r abg 2 CQt / x coth ( 6s ) _ ! _ ( 15 ) 

7^ + M Jo (Airsy 

Here s is the Schwinger proper time, tr {T a T a ) = Cf • trljv c xiV e = — l)/2, and Cf = 
(TV 2 — 1)/(2N C ) is the quadratic Casimir operator of the fundamental representation of the 
SU(iV c )-group. Furthermore, a and 6 can be expressed through the electric and magnetic 
fields corresponding to the Abelian strength tensor as 

1 



a 2 



E 2 — H 2 + v / (E 2 -H 2 ) 2 + 4(E-H) 2 
H 2 - E 2 + v/(E 2 -H 2 ) 2 + 4(E-H) 2 



2 

The effective action represents the sum of infinitely many diagrams with one quark loop and 
various numbers of external .A^-lines. Due to the factor e~ M2s , one can readily see that the 
parameter of the expansion of the effective action in the number of external lines is (G 2 ) /M 4 . 
Therefore, in the heavy-quark limit under study, the leading contribution is produced by 
the diagram with only two external lines. One can further bring together Eqs. f lT0|) - f|T5|) 
and proceed, when integrating over the J-^-components in Eq. (TlT!) . to the 6-dimensional 
spherical coordinates. Using the value of the 6-dimensional solid angle, Q e = tt 3 , one arrives 
at 

4n c n, i .^r,*^ (16) 



heav y (8ttC) 3 (4tt) 2 3 



o 



(g 2 ) 

where C = and B is the radial coordinate in the 6-dimensional space. In the heavy- 

quark limit at issue, the factor e~ M2s filters out all proper times but those for which 
s 1 ~ M 2 . For this reason, the upper limit in the 5-integral can be replaced by +oo, 
thereby decoupling the s- and the 5-integrals from each other. This yields the heavy-quark 
condensate 

2N C N { C 

Wheavy = ^ 
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which reproduces correctly the known result of QCD sum rules By calculating the 
-B-integral exactly as 

jf 1 " dBB 7 e~^ = 3(8C) 4 - ^ (l + 2ACs 2 [l + lQCs 2 (l + 8Cs 2 )] \ exp (-l/(8Cs 2 )) , 



one can see that corrections to Eq. ffTTj) . arising from the entanglement of the two integrals 
in Eq. (Tl6|) . are suppressed by an exponentially small factor 0(e~ sc ). 

Let us now proceed towards the limit of small current quark masses. To this end, we 
model excitations of the minimal area by extending the gauge field as 

A%^Al = A%- -a%. (18) 

Let us define the propagator (g 2 ~o'p(x)~o'p(y)) as K^ u (x — y). Then, this correction 
will be equivalent, in the lowest-loop approximation, to an extension of the action 
fx [^(Tm-^m + m )"0] by the current x current interaction through the kernel K /Jtl/ (x): 

S= /[^(7^ + m)V]+ / j;(x)K, v (x-y)jZ(y). (19) 

J x J x,y 

Here = — igT a A^ is the covariant derivative with A a ^ yielding the minimal-area law. 
We furthermore assume K^x) to be long-range. Then, the current x current interaction 
mediated by the long-range kernel K^ v {x) should be thought as the generator of longitudinal 
string vibrations, which are not accounted for by the minimal-area law alone. Our aim is 
to demonstrate how a long-range interaction kernel K^x) can be absorbed into the scale- 
dependence of the constituent quark mass m, as it was the case for the GNJL model - see 
Eq. (j2J). We note once again that, for heavy quarks whose trajectories fit into the circle 
of diameter a, the constituent quark mass goes over to the current mass M. Instead, light 
quarks acquire their constituent mass primarily in the course of their motion at distances > a. 
This observation leads us to the natural assumption that constituent quark masses should 
not depend on a local point at the quark trajectory, but rather on the global characteristic of 
the trajectory - in this case, its mean size d. This assumption is translated to the world-line 
formalism by means of the proper-time dependence of the constituent mass: m = m(s). 

To see this, let us use the known world-line representation of the quartic self-interaction 
of a matter field. In the simplest case, when the role of j" is played by the scalar density of 
a complex-valued Abelian Higgs field, this representation, reviewed in Appendix A, involves 
an auxiliary scalar field a(x). This field is used to disentangle the quartic self-interaction 
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via the Hubbard-Stratonovich transformation. In the actual case under study, with the 
non-Abelian vector current reinstated, this a-field acquires a Lorentz and a color indices, 
that is becomes <r". A counterpart of Eq. ( 1A3I) . J V^jjVijje~ s ~ e~ r ^, defines an effective 
action rL4"]. Its representation analogous to Eq. (1A4I) reads: 

r[A% = J VaZe-y^-^nA",*^ (20) 

where T[A°, <r£] = -tr In [y^D^ + iT a o") + m] . Comparing T[A a ^ a«] with T[A^, a] , we see 
the main difference of the fermionic case at issue from the bosonic one of the Abelian Higgs 
model: Because in the fermionic case the kinetic term is first-order, the field a?, is included 
directly to D^, and not to D 2 as in the bosonic case. Notably, in the fermionic case, the 
field cr" just extends the gauge field A ^ according to Eq. ( TT8|) . In terms of such an extended 
gauge field and its strength tensor J 7 ^, = d ^A^ — d u A a ^ + g f abc A^A^, a fermionic counterpart 
of Eq. (1A5I) reads 

T[A«,a;} = -2N { V f° -e""^ / Vz, [ V^e~ & W [ Zfi , ^, A%. (21) 

JO s J P J A 



Here 



W[z„,%,A%] =trV exp 



f drT a {Alz,-^^ v T; u ) 
Jo 



(22) 



V is the 4D volume occupied by the system, P and A stand, respectively, for the periodic 
and the antiperiodic boundary conditions f p = J z ^ =z ^ and J A = ^ = _^ r y an d 
the trajectories z^{j) obey the equation f° drz^ij) = 0, which constrains the center of a 
trajectory to the origin. Furthermore, the same lowest-loop approximation, used to derive 
Eq. (1A6I) . now yields 

J Va^e-y*^' 1 ^^ Wlz^^A^] ^ Wlz^^A^] ■ e ~c t f d Zfl f dzlK^-*') . (23) 

Thus, similarly to the bosonic case, Eq. (1A6jl . in the fermionic case we have also arrived at 
the self- interaction of the loop mediated by the kernel K liv {x). 

We are now ready to show that a long-range kernel can effectively be substituted by some 
proper-time dependent constituent quark mass m(s). That is, we are seeking such a function 
m(s) as to obey the equality 

J VZffi~^^ dTzl-m 2 (s)s _ J J) z ^ e -\ Jo drzl-Ct § d Zf , § dz' v K„ u (z-z') _ ^4) 
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By way of illustration, let us consider the simplest, quadratic, kernel K^ v (x) = -^-5^ v x 2 ^ 
where K is some parameter of the dimensionality of mass. Since for such a kernel one has 
§ dz^§ dz' v K^ u {z — z') = —Kq (</> dz v Zy) A \ the world-line integral on the right-hand side of 
Eq. ( |24"1) can again be calculated by introducing some auxiliary constant antisymmetric field 
B^ v as follows: 

-C f § dzp § dz^K^z-z') _ TT / dB I p~W~5 B ^ So 'drz^zv 

where for the rest of this section 

C = QK*. (25) 

The resulting world-line integral is given by the bosonic Euler-Heisenberg-Schwinger La- 
grangian: 

f Vz e ~ IS = — . (26) 

J P ^ (47rs) 2 sin(6s) sinh(as) 

Expanding this Lagrangian up to the leading nontrivial term as 

abs2 ^1 s2 y B \ 

sin(6s) sinh(as) 6 ^ v 



and proceeding in the B^-integration to the 6-dimensional spherical coordinates, we obtain 

1 -m\s)s 1 1 f 1/i8CS2) . 2 A \ 

e m (s)s ~ — — - ■ - / dyy 2 1 — y e y , 



(4vrs) 2 (4tts) 2 2 J 



where y = (8C) 1 ^ -BjL is the radial integration variable in the 6-dimensional space. 

~ ,2/ 



Finally, introducing instead of s a dimensionless proper time A = vCs, we obtain m (A) in 
the units of \/C: 

m 2 (A) I 

•in , 

2./o V 3 



-- - In 



1 /- 1/(8A \ 2 A 4A 2 , _ (/ 
dyy 1 — r-y e y 



/(A). (27) 



^ A 

This function is plotted in Fig.[TJ It has a maximum at A extr ~ 0.62, above which m 2 steadily 
falls off and can be approximated by the function 



Thus, at sufficiently large proper times relevant for the spontaneous breaking of chiral sym- 
metry, the effective constituent quark mass squared approximately falls off as 



13 



12 




10 - 



2 



8 



6 



4 












5 



10 



15 



20 



FIG. 1: The dimensionless constituent quark mass m 2 j\[C as a function of the dimensionless 
proper time A = \^Cs, Eq. (|27p . The maximum of m 2 /y/C is reached at A cx tr — 0.62. At A > A ex t r , 
is also shown the approximating function Eq. (|28p . 

We can now obtain the effective constituent mass as a function of the physical mean size 

n 

d of a trajectory. It has been shown in Ref. [3J that the classical trajectories of a confined 
light quark are circles, so that the mean size d is the diameter of a circle. We notice that, 
for an s-dependent mass m(s) at issue, the length of the quark trajectory L is represented 
by the integral L(s) = m(r)dr. Using the function /(A) introduced in Eq. fl27j) . we can 
write the length L as an integral over the dimensionless proper time A as 



We can furthermore use the result of Ref. [3] that a light-quark trajectory has Hausdorff 
dimension 4, that is its length grows as the fourth power of its mean size, L oc d 4 . Since 
the minimal mean size of a trajectory, for which the quark is still confined, is d m i n = 2a, we 
obtain the following formula for the length of a trajectory in terms of its mean size: 



Equating to each other the two expressions for the length, Eqs. (1291) and ( )30l) . we can obtain 
the dimensionless proper time in terms of the mean size of the loop, X(d). We choose 
the range of distances d G [d min ,2fm], use the commonly accepted value of the vacuum 
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FIG. 2: The dimensionless proper time A = \/Cs as a function of the mean size d of a trajectory. 



correlation length 19h21| a = 0.34 fm, and evaluate the parameter C, Eq. ( 1251) . with the 



strength of the kernel Kq = 300 Me V (cf. the second reference in [6|]). We notice further 
that X(d) is a monotonically increasing function, and that already A(<i m i n ) ~ 0.79 exceeds 
the position of the maximum in Fig. [H A ex tr — 0.62. This fact facilitates the procedure 
of the numerical solution of the equation L(X) = L(d). Namely, instead of using in the 
length-functional, Eq. ( |29|) . an integral representation of the function /, Eq. ( |27l) . we use 
the approximating function / appr , Eq. (1281) . In this way, we obtain the function X(d) and 
plot it in Fig. [2j Substituting now A in Eq. ( 1271) by this function A(<f), we obtain also 
the effective constituent quark mass in terms of the mean size of a trajectory. Such a 
function m(d) is plotted in Fig. [31 In particular, we find that m(d) acquires its commonly 
accepted phenomenological value of 300 MeV at d ~ 1.30 fm. The maximum value, which 
the constituent mass acquires, is the one at d = d m i Q , that is m max = 953 MeV. Thus, we 
have shown that, for a given confining kernel K^ v (x), we also have, like in the GNJL case, 
a running mass m(d), this time a function of the loop size d. In particular, the constituent 
quark mass m max defines the scale for the highest possible replica, arising when the loop is 
shrunk down to the size of d m i n . We notice finally that, in addition to d min , there always 
exists some limiting maximum value <i max . It is associated with the string-breaking, as at 
d > d max the system becomes unstable due to the quark-antiquark pair creation. Lattice data 
and analytic models suggest some value around 2 fm 2^] for the string-breaking distance. 
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FIG. 3: The effective constituent quark mass m as a function of the mean size d of a trajectory. 

We therefore take the combined theoretical evidence for the existence of replicas, coming 
from both the GNJL and the effective-action formalisms, as a hint that they might be 
realized in Nature. However, in order for this to have a chance to happen, we still have 
to evaluate, through the use of the Tolman-Oppenheimer-Volkoff equation, the stability 
condition for a replica-filled domain versus the gravity field it might produce. To do this, we 
need to know the equation of state for such macroscopic objects. This will be the subject of 
the next section. Finally, in section IV, we extend the effective-action approach to evaluate 
the leading gravitational correction to the quark condensate. There, it will be shown that for 
distances larger than the hadronic scale, gravitational corrections to the flat-space results are 
negligible. This fact allows us to use the standard flat-space value of the chiral condensate 
to study the equation of state of such macroscopic domains. 

III. THE ENERGY DENSITY AND MACROSCOPIC PROPERTIES OF A 

REPLICA-FILLED DOMAIN 

Let us proceed with the evaluation of the full internal-energy density e(T) of a macroscopic 
domain whose quantum states are built on top of a replica state. The full internal-energy 
density of such an object is given by the sum of the replica-vacuum contribution and the 
contribution produced by hadronic excitations above that vacuum (clearly, staying in the 
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lowest-energy domain, we restrict ourselves to pionic excitations): 



e(T) = £ yac (T)+ £h (T). 



(31) 



We consider first the vacuum contribution e vac (T) . It is defined as a difference of the internal- 
energy density of the excited vacuum, corresponding to a replica state inside the domain, 
and the internal-energy density of the unexcited vacuum outside the domain: 



evacCO = e R {T) - e (T). 



The internal-energy density of the unexcited vacuum has the form 



(32) 



eo(T) = - 



32tp 



(G 2 ) T + (m u + m d ) 



(33) 



where the one-loop /3-function coefficient for N c = 3, N{ = 2 is b = 29/3, and the ap- 
proximation [23] (ijrij)} ~ (uu) ~ (d(£j has been adopted. One can further use the known 
temperature-dependent gluonic and chiral condensates, which at temperatures T <C m n of 



interest read 



23) 



g 2 ) t = <g 2 > - (^) , <^> T = (m 



3m^T /m, 

J 7T 



(34) 



Here, (G 2 ) and (ipvp) are the zero-temperature values of the condensates in the unexcited 
vacuum, f n is the pion decay constant in the same vacuum, and we denote the sums entering 
virial expansions of the grand canonical ensembles of relativistic bosons as 

KJnx) 



Su(x) = 



n=l 



n" 



where K v is the MacDonald function. Inserting Eqs. ( 1341) into Eq. (1331) . and using the 
Gell-Mann-Oakes-Renner relation |8| (f^m^) 2 = —{m u + m^) (ipip), we obtain 



Eo( r) = £o (o) + ^ Sl (^ 



(35) 



where Eq(0) is the internal-energy density of the unexcited vacuum at zero temperature. For 
the internal-energy density of the excited vacuum, in which the Gell-Mann-Oakes-Renner 
relation holds equally well (cf. the last reference in 13J), an expression similar to Eq. ( 135]) 
reads: 



e R (T) = e R (0) + 



3m! T 



8tt 2 



Si 



m 



(36) 
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where m 7TR is the mass of the pion in the replica. Equations (152]) . (|35|) . and ( 136]) yield 

^ac(T) = e + g [m 3 ^ (^) - mlS, (!£)] , (37) 

where e = £#(0) — £o(0). 

We find now the hadronic contribution £h{T) to the internal-energy density (13T]) . By 
neglecting the outer temperature with respect to the temperature inside the domain, we 



neglect a 



reads 



23| : 



so the outer pressure. The pressure of the relativistic pionic gas inside the domain 



P^H^&H?). (38) 



2tt 2 z V T 

The hadronic contribution to the internal-energy density can be found by the standard 
thermodynamics formula 

e h (T) = T^p--p h (T). (39) 

Finally, according to Eq. ( |3T]) . the full internal-energy density is given by the sum of Eqs. ( l3Tj) 
and ( 139]) . Noticing that m nR > [3], we see that the leading correction to the formula 
e(T) ~ e at temperatures T reads 

Higher corrections are getting more important with the increase of temperature. Using the 
values [kJ 

m^ = 250 MeV, (^)fi = -(100 MeV) 3 , e~(250MeV) 4 , (41) 

and the standard 7r-meson mass m n = 140 MeV, we evaluate numerically the internal-energy 
density (|3T]) along with its approximate expression ( 140]) . and plot these quantities in Fig. HI 
We observe, in particular, the complete numerical irrelevance of the hadronic contribution 
up to the temperatures ~ 25 MeV. For this reason, in what follows we disregard the hadronic 
contribution altogether, and approximate the internal-energy density of the domain by the 
constant value 

e(T) ~ e = (250 MeV) 4 . (42) 

It follows from Eq. ([6]) that, with an exponential accuracy, a replica-filled domain oc- 
cupying a sufficiently large volume is stable against the decays into the ground-state 
vacuum. The phenomenon providing this stability is the very strong correlations between 
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FIG. 4: The full internal-energy density, e(T) • 10 9 (MeV) 4 , given by Eq. (|3ip . and its ap- 
proximation, Eq. (|40p . The constant value, s(T) ■ 10~ 9 = 3.90625 (MeV) 4 , corresponds to 
e(T) ~ e = (250 MeV) 4 , Eq. ([42]) . 



pions in the replica. We notice also that, with the energy density (j4"2j) . the mass of a 
replica-filled domain grows with the volume as e x V^ 3 \ so that the problem of stability 
with respect to the gravitational collapse becomes relevant. For a constant energy- density, a 



solution to the so-called Tolman-Oppenheimer-Vo 



koff equation, which provides the metric 



of a spherically-symmetric object is well known 12j. In Appendix B, we briefly discuss this 
solution, which defines the maximal possible radius of the domain [see Eq. (IB 8ft ] as 

1 



i? G =^==^14km, (43) 
V oireG 

where G = 6.7- 10 -39 GeV -2 is the gravitational constant. Notice that the obtained estimate 
for the maximum radius is of order of the typical radius of a neutron star, that is about 
10 km. This analogy extends further to the rotating neutron stars (pulsars). Namely, 
consider a domain of mass M rotating with the angular velocity Q. Then, for an element 
of matter of mass \x situated at the equator of the domain (where the centrifugal force is 
maximal), the Newton's second law has the form 



where the correcting factor £ ~ 0.65 accounts for the effects of General Relativity 24]. 
Noticing that M = AtceRq/S, we obtain from this Newton's equation the period of the 
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domain's rotation 

T G = — ~ 0.6 msec. (44) 

Notice that periods of rotation of certain pulsars can be as small as milliseconds, that 
qualitatively agrees with the value (1441 . 

We end this section with a brief discussion on the stability of the replica-filled domains 
against the decays of their pionic excitations into photons. Such anomalous decays appear 
suppressed compared to the similar conventional decays in the unexcited vacuum. To see this 
suppression, we use the Gell-Mann-Oakes-Renner relation jg], which yields the pion decay 
constant f n . Since in the replica the chiral symmetry is "less broken" than in the unexcited 
vacuum, the absolute value of the chiral condensate in the replica is smaller, whereas the 
pion mass is larger. Using the values of m nR and (i/jiP)r quoted in Eq. (141jl . one can readily 
evaluate a two-photon decay width of a "replicated" pion n R as 

T(n R 77) = r(vr 77) " ^ = 77) • { JT^ ^ * 0.17eV, (45) 

whereas T(7r — > 77) ~ 8eV. The smallness of the calculated decay width r(7i\R — > 77) 
suggests that a replica-filled domain would appear to an external observer as a dark massive 
body made of strongly correlated pions. 



IV. A GRAVITATIONAL CORRECTION TO THE CHIRAL CONDENSATE 

The estimate of the radius of a replica-filled domain, made in the previous section, em- 
ployed the flat-space result for the vacuum energy-density. It is clear, however, that chiral 
symmetry breaking experiences influence of gravity and that, in general, the MGE and the 
Tolman-Oppenheimer-Volkoff equation are mutually coupled. On the other hand, the esti- 
mated radius Rq of a stable domain constitutes a few kilometers, while the scale of chiral 
symmetry r x is as small as just a few fermi, that is many orders of magnitude smaller than 
Rg- In other words, we have a well pronounced hierarchy of scales in the problem, and thus 
we anticipate that gravitational corrections to chiral symmetry are negligibly small. 

In this section, we justify the above qualitative argument by evaluating the leading cor- 
rection acquired by the chiral condensate in the gravitational field of a replica-filled domain. 
Such a correction is expected to be of the form 0(d 2 )/m 2 , where m is the constituent quark 
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mass, and the (9(d 2 )-terms contain in particular the Ricci scalar TZ. Within the effective- 
action approach, the quark condensate stems from the formula (ipil>) = — (T[Afy. In the 
flat space, yet another derivation of Eq. (fT6l) can be performed (cf. Appendix A of Ref. Q]), 
which is based on the one-loop effective action 



Vz„ 



wher e ag ain C 
Ref. 



Vijjf, exp 



2NM 



- / dr 



ds 
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\fl<U 



+00 



dB 



e~w x 



1 1 • i 



(Atts) 



, 1(iV , and B^p is an auxiliary constant Abelian field. As is shown in 
the result of the functional integration in this formula is given by the Euler- 



Heisenberg-Schwinger Lagrangian, which can further be Taylor-expanded in s: 

1 



{••■} 



(Atts) 



(46) 



n=0 



That is, of course, just the heat-kernel expansion of the (logarithm of the) determinant of the 
Dirac operator of an electron in the electromagnetic field B^ v . Each term of this expansion 
is described by the one-loop diagram with the corresponding number of external lines of 
the -B^-field. In particular, the leading non- vanishing term, with the coefficient = \Bj lu , 
which corresponds to the diagram with two such lines, recovers Eq. (TIB"]) , where B = ^ B 



[J.V ' 



A distinguishing property of the heat-kernel expansion in the curved (Euclidean) space 
is that some of the diagrams with two external lines of the -B^-field contribute not only to 
the 2-nd, but also to the 3-rd term in Eq. (146 p . That is because such diagrams additionally 
describe an interaction of the electron with the gravitational field. This interaction involves 
operators of dimensionality (mass) 2 , such as the Ricci scalar TZ, which are accompanied by 
an extra power of s. Explicitly, the coefficient 03 in the presence of the gravitational field 



has the form 
1 



25] 



03 



360 



[bTZB^B^ - 2&W"'B tlX B v A + 2U^ Xp B^B Xf> + 2A{V v B^){V x B tlX )] . (47) 



Notice also that the measure of integration over the B^-field is changed in the curved space, 
where 
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Performing the s-integrations, we arrive then at the following intermediate expression for 
the quark condensate: 



N c N f (det g^ v fl 2 



471-2 



n 

\fl<U 



dB 



\ PbX p (a 2 a 3 



- + —,h (48) 



where a<i = ^B^B^, while 03 is given by Eq. f H7|) . In particular, for 03 = and g^ v = 5^, 
Eq. (TJH]) reproduces Eq. f lT7|) . The .B^-integration in Eq. f )48|) can be performed through 
the formula 



(det# 



fiv\3/2 



(8ttCY 



n 



dB 



B x PB X 



hp 



We furthermore proceed to the Minkowski space, since the metric of the domain to be 
eventually used is Minkowskian. There, due to the constancy and antisymmetry of B^, one 
has V V B^ LV = B^ v d v In y/—g, where g = detg^. Finally, contracting the indices, we obtain 
for the quark condensate the following result: 

1 



[W) = <^> 



(49) 



where \ip^/ is the condensate in the flat space. 

Equation ( H9|) can now be used to evaluate the quark condensate in the gravitational field 
of a replica-filled spherical domain of a constant energy- density. We start with the inner part 
of the domain and assume that it has the maximum possible radius (|43|) . Now, inserting 
Eq. (B7) into Eq. (B9), we obtain the following expression for the temporal component of 
the metric inside the domain: 

8/9 



9oo{r) = - • exp 



dz' 
l-z' 



1 + 



3^1^ - 1 
1 - VT^z 7 



(50) 



where z = 8ireGr 2 /3 is such that z < 8/9. Using also the other components of the metric, 
namely 

' -r 2 sin 2 #, (51) 



9u 



z-1 



922 



-r , £33 



we have for the corresponding term in Eq. ( l4~9l) : 



g^dAnV^g-dAn- 



(z-1) 



cot 2 fl 



(52) 
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The differentiation in this formula is straightforward and yields 



g^dAn^-d u \n^ 



z - 1 



- cot 2 e 



(53) 



Due to the spherical symmetry of the domain, we are interested only in the radial dependence 
of the gravitational correction to the quark condensate, while averaging out its angular 
dependence. That can be done by noticing that the azimuthal angle enters the invariant 
volume element as d^x^—g oc <i#sin 2 #, thereby defining the following normalized angular 
average: (• • -) e = - J* dOsm 2 9(- • •). In particular, the so-averaged angular-dependent term 
in Eq. fl53|) yields just unity: (cot 2 ^)^ = 1. 



The Ricci scalar 1Z entering Eq. 



and b(r) = ln(— gn) as follows [12 
K = e 
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(54) 



A straightforward calculation yields for it a simple result: 

rr, 12(1-^) 

n= K r2 J . 

The full 6*-averaged correction to the quark condensate, Eq. (JUD, thus reads 



(55) 
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The obtained correction to the absolute value of the quark condensate, | 
definite, and it gets smaller than 1 when the following condition is met: 



- - 2z 

J (56) 
is positive- 



in 

z 



+ 



2 < 



15m 



, , _ (57) 

To numerically evaluate the obtained correction in the chiral limit, we use the internal- 
energy density ( T42l) . as well as the values of the constituent quark mass m = 300 MeV 
and the gravitational constant G = 6.7 ■ 10~ 39 GeV -2 . We obtain for the right-hand side 
of inequality ( 157)) : 15m 2 /(47reG') ~ 4.1 ■ 10 39 . Seeking distances r at which the obtained 
correction to the condensate can become relevant, we get from the definition of z and Eq. (J5j 



r < 



0.55 fm. 



m 



This entails the following correction to the quark condensate: 

<^> = (H>)o [l + 0((r x /R G ) 2 )} , 



(58) 



(59) 
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where Rq is given by Eq. (143]) and r x = 1/m, with m being the effective constituent quark 
mass, which appears due to chiral symmetry breaking. One can thus see that, as was 
anticipated, inside the domain, the effect of its gravitational field on the chiral condensate 
is negligibly small. 

Outside the domain (r > Rq), the metric is Schwarzschild (cf. Appendix B), therefore 
1Z = 0. The remaining term in Eq. ([49]) reads for this metric: 



g^dAn^-d v \n- 



1(^-1 

r 



cot 2 d 



(60) 



where r g = 2GM is the Schwarzschild radius, with M being the full energy of the domain 
- see Appendix B. Once averaged over 6, Eq. (160]) yields the following positive-definite 
correction to | (i>i)) \ - 



(U-Q- (g^ In ^ ■ d v In yf^ > ) 



4(l-^)+l 

r 



(61) 



60m 2 v w a 1 v * v v /eJ 10 ( mr )2 

Noticing that r g /Rc = 8/9, we see that [4 ( 1 — + l] grows from the value 13/9 at r = i?c 
to the asymptotic value 5 at r — > oo. Therefore, we can estimate the obtained correction as 
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1.4-10 



-39 



(62) 



Such a numerical smallness of the result originates, of course, from the value of the gravita- 
tional constant G. We conclude that, outside the domain, the effect of its gravitational field 
on the chiral condensate is negligible as well. 



V. SUMMARY 



In this paper, we considered the possibility of formation in the early Universe of macro- 
scopic domains filled with coherently excited 7r-mesons. At temperatures below 25 MeV, the 
energy density inside the domains is primarily defined by the energy gap over the physical 
QCD vacuum of the first excited coherent state (the so-called first replica) in the GNJL 
model. Instead, the pionic excitations above the replica are shown to be suppressed at such 
temperatures. In this way, using the constant energy- density ( 1421) . we obtained the maxi- 
mum possible radius of the domain, Eq. ([431 . which turns out to be similar to the radius of 
a neutron star. Furthermore, for rotating domains, the period of their rotation (T44|) is also 
similar to that of some class of pulsars. 
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Since our approach is based on the energy density in the flat space, we have checked 
its self-consistency by demonstrating a smallness of the leading correction to the flat-space 
chiral condensate appearing in the gravitational field of the domain. To this end, we have 
first illustrated the correspondence between the GNJL and the effective-action descriptions 
of spontaneous chiral symmetry breaking. Then, using the effective-action approach in 
the curved space, we derived the leading (that is, containing the minimal number of the 
derivatives) gravitational correction to the quark condensate, Eq. fj4*9|) . Finally, using the 
gravitational metric of a const ant-energy- density domain, we evaluated the correction to the 
chiral condensate both inside and outside the domain [cf. Eqs. ( 1561) and (16T1)]. Outside the 
domain, this correction is extremely suppressed due to the smallness of the gravitational 
constant G [cf. Eq. ( |62l) ]. whereas inside the domain the correction could only get relevant 
at the distances to the center of the domain as small as those given by Eq. fl58l) . Such 
distances are, however, 20 orders of magnitude smaller than the macroscopically large radius 
of the domain, which makes the obtained correction numerically irrelevant. The smallness 
of the correction, that the chiral condensate receives in the gravitational field of a domain, 
completely validates our use of the flat-space chiral condensate. 

In conclusion, we put forward a conjecture that the domains of coherently excited pions 
could have been created in the early Universe. Such domains are hot (as the estimates 
made in this paper remain valid up to temperatures of order 25 MeV) and stable against 
the gravitational collapse up to the maximum radius of about 14 km. Moreover, since the 
decay width of the coherent pionic states into photons is negligibly small [cf. Eq. (1451)]. 
these domains cannot evaporate by means of the electromagnetic radiation, which leads to 
their thermal insulation from the outer world. Since one can also argue for the stability of 
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14] . such encapsulated 



the coherent pionic states against the strong and weak decays 
domains can have had a chance to survive till the present time, remaining however dark to 
external observers. 
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Appendix A: A world-line representation of the quartic self-interaction of matter 

fields 

In this Appendix, we recall the known world-line representation of the Abelian Higgs 
model with a non-local quartic interaction. Such a model is a bosonic prototype of the 
GNJL model. The part of its action containing the complex-valued Higgs field <&(x) reads 

S= [ (\D^\ 2 + ™ 2 m 2 )+ [ \*(x)\ 2 K{x-y)\*(y)\ 2 , (Al) 

J x J x,y 

where = <9 M — igA^ is the covariant derivative. To integrate over the Higgs field, one ap- 
plies to the quartic term the Hubbard-Stratonovich transformation. Since it is the absolute 
value |$|, which enters the quartic term, the auxiliary scalar field cr(x) is real- valued: 



e -L, y \Hx)\ 2 K(x-y)\<S>(y)\ 



J Vae-Ui^- 1 *^ 2 ), (A2) 



where K~ x is a local operator, whose Green function is K, that is K~ l (x)K(x) = 5(x). (In 
particular, for an ordinary Abelian Higgs model with the local quartic interaction K(x) = 
k5(x), the operator K~ l is just a number: K~ l = 1/k.) Integration over the Higgs field 
then yields 

J T?$D$*e~ 5 = J Vae-^^ aK ' la - r[A ^ a] ~ exp 1^- J Do e - ^ L aK ^ Xrj TfA^, a] j- , (A3) 

where r[A M ,cr] = tr ln(— D 2 ^ + m 2 + io) is the effective action. The last approximation in 
Eq. flA3|) was made analogously to the one-loop approximation for the effective action with 
respect to the gauge field (A^ in the present example), adopted in the bulk of the paper. It 
means that the effective action describes a loop of the Higgs field with an arbitrary number of 
external lines of the gauge field, but does not describe gauge-field exchanges inside the loop 
and/or gauge-field interactions of two and more such loops. Similarly, Eq. (IA3j) assumes the 
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lowest-loop approximation for the effective action with respect to the cx-field. We consider 
now the averaged effective action 

= J Vae-y^^TlA^a] (A4) 
and apply the world-line representation to write [cf. Eq. ( 121 J) ] 

T[A^a] = 2V n^ e -™*s f Vz ^-f°dT[±zl+i<j(z(T))] W [ Z ^ A ^ (A5) 
Jo s J P 

Here, W[z^,A^\ = exp (ig J* Q S drA^z^j is the Abelian phase factor along the loop, and the 
overall factor of 2 is because of the two (real and imaginary) components of the Higgs field. 
Integrating now over the cx-field, one obtains the following effective action: 

T[A^]=2V r^e-" 1 ' 3 [ Vz^-yo^l-I a dr^dr'K[ Z (r)-z(r')) W ^ A ^ (A6) 

Jo s Jp 

As one can see, because of the two-point interaction produced by the nonlocal kernel K(x), 
we have arrived at the two-loop effective action (instead of the one-loop effective action, 
which one gets in the same lowest-loop approximation when accounting only for the interac- 
tion of the Higgs field with the gauge field), which cannot be calculated analytically In the 
particular case of a short-range (or a local) kernel, the term J Q S dr J Q S dr'K[z(r) — z(t')\ sup- 
presses in the world-line integral self-intersecting trajectories. For this reason, the Abelian 



Higgs model with a short-range nonlocal quartic interaction is known to 
representation for the ensemble of self-avoiding bosonic random walks 



3e a field-theoretical 



27]. 



Appendix B: The gravitational field of a spherically-symmetric object with a 

constant energy density 

In this Appendix, we summarize some known facts [12| about the gravitational field, 
the energy density s, and the pressure p of spherically-symmetric objects, with a par- 
ticular emphasis on those with the constant energy density. Throughout this Appendix, 
for terminological simplicity, we will call such objects stars. Usually, one assumes that 
the matter forming a star is a perfect fluid. This implies the energy-momentum tensor of 
the form T^ v = (p + e)u^u u — pg^ u , with u^{x) being the four-velocity of the fluid, such 
that g^yU^u 11 = 1. In the local rest frame of the fluid, where = ( v /<?oo, 0), the energy- 
momentum tensor is diagonal: 

v = (p + e)u^u v - pd** v = diag (e, -p, -p, -p). 
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Furthermore, since the star is spherically symmetric, the x-dependence of e and p is reduced 
to their dependence on the spatial distance to the center, r. For this reason, one uses the 
3D spherical coordinates (r, 9, (f>) to find the metric components (700 = e a ^ and gw = —e b v\ 
while the two other components read #22 = —f 2 and #33 = — r 2 sin 2 9. Outside the star, that 
is for r > R, where R is the star radius, the energy-momentum tensor T^ v vanishes, and 
one obtains the Schwarzschild metric 

a (r) = -b(r) = In (l - ^) (Bl) 

with the Schwarzschild radius r g = 2GM, where G is the gravitational constant, and M = 
An $q drr 2 e(r) is the full energy of the star. 

For r < R, it is convenient to introduce the energy inside the radius r, given by Ai(r) = 
4vr ^dr'r' 2 s{r'). The Einstein equation TZ° - \U = 8nGT° reads 

e- b fdb 1\ 1 o ^ 
r \dr r J r z 

Its solution b(r) = — ln(l — ) clearly goes over to the function ( 1B1I) at r = R. The 
function a(r) inside the star can be found from the covariant conservation of the energy- 
momentum tensor, V ^T^ u = 0, which explicitly reads 

- d,p ■ gT + 0„ [(p + e)«"t*1 + (p + e) (rlu x u» + I^tiV) = 0. (B2) 

One assumes furthermore the hydrostatic-equilibrium condition, which implies the Xo- 
independence of not only p and e but also of = (-y/poo, 0). In particular, this condition 
yields d^p + e)u^u u ] = d [(p + e)u°u°] = 0. Multiplying the rest of Eq. ( 1B2I) by g pu , one 
has 

^ = rW + ^ttV. 

The first term on the right-hand side of this equation vanishes, while the second term leads 
to the desired equation 

d\ng m = _ 2 _ dp/dr 
dr p + e 

The solution to this equation, 



9oo(,r) = g 00 (R) ■ exp 



R 



dr 



, dp/dr' 



with 0oo(i*) = l-^| (B4) 



p + e 

defines the temporal component of the metric inside the star, as well as the function a(r) = 
\ng 00 (r). Similar to the function b(r), such a function a(r) goes over to Eq. fIBlj) at r = R. 
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The equation for p(r) can be obtained from the Einstein equation 1Z 1 1 — ~1Z = 81TGT 1 1( 
which is equivalent to 



da 



1 



_ = e " ( - + 8nGrp 



1 

r 



To this end, it suffices to use da/dr given by Eq. (IB3I) . 4^ 



, dp / dr 



and the relation 



1 — 2GM./r) 1 . That yields the so-called Tolman-Oppenheimer-Volkoff equation for 



p(r) [Tj 



dp 
dr 



GeM 



1 + 



Airr 3 p 
M 



(B5) 



Together with the equation ^ = 4?rr e and the equation of state, it forms a set of three 
equations for the three unknown functions, that are p, e, and hA. Furthermore, plugging 
Eq. flB5j) into Eq. (IB4j) . one also obtains the temporal component of the metric inside the 
star in terms of the functions p and Ai: 



9oo(r) = 0oo(-R)exp 



-2G 



R dr' 



-f'2 



20M\~ 1 
r' J 



[M + ATcr' 3 p) 



(B6) 



In the particular case e = const, one has M = 4fer 3 , and Eq. ([B5} can be integrated 



analytically to yield 



p(r) 



1 - f eGr 2 



1 - f eGR 2 



3J1 - ^-eGR 2 



1 - f eGr 2 



(B7) 



where the boundary condition p(R) = has been taken into account. For the denominator in 
this formula not to vanish at any r < R, one imposes the requirement 3a /l — ^-eGR 2 — 1 > 0, 



which defines the upper limit for the star radius: 

R< 1 



(B8) 



v3neG 

Finally, according to Eq. (1B6|) . the temporal component of the metric inside the star of a 
constant energy density reads 

•R dr > r i 



9oo{r) 



1 



R 



exp 



-87TG 



j _ 8TTsG_ r /2 \3 



+ p(r'] 



(B9) 



where the function p(r') is given by Eq. (1B7j) . 
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